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Abstract;  A  robust  procedure  for  location  and/or  scale  parameter 
estimation  is  presented  which  utilizes  the  asymptotically  best 
linear  unbiased  estimators  (ABLUE's)  based  on  k(*N)  of  the  N  sample 
quantiles.  Using  regression  design  techniques  a  method  is  developed 
for  selecting  sample  quantiles  which  furnishes  the  corresponding 
parameter  estimates  with  good  robustness  properties  relative  to  a 
given  finite  set  of  known  probability  laws.  The  problems  of  robust 
quantile  selection  for  the  estimation  of  a  particular  population 
quantile  or  in  the  presence  of  left  and/or  right  hand  censoring  are 
also  considered.  - 

1.  Introduction. 

A  location  and  scale  parameter  model  assumes  that  a  random 

sample,  X^,...,^,  has  been  obtained  from  a  distribution  of  the  form 
X-u 

F(— — )  where  u  and  a  are  respectively  location  and  scale  parameters. 

An  important  problem  associated  with  this  model  is  the  development  of 
parameter  estimates  which  behave  robustly  (e.g.,  insure  high  asymptotic 
relative  efficiencies)  over  several  or,  perhaps,  many  choices  for  the 
form  of  F.  In  this  paper  we  provide  an  approach  for  obtaining  such 
estimators  which  utilizes  the  asymptotically  best  linear  unbiased 


estimators  (ABLUE's)  of  y  and  a  based  on  k<N  sample  quantiles. 


The  ABLUE's  are  easily  computed  estimators  whose  properties 
have  been  extensively  studied  (c.f.  Ogawa  (1951),  Sarhan  and 
Greenberg  (1962),  Eubank  (1981)  and  the  references  therein).  Let 
Q  denote  the  sample  quantile  function  defined  by 

Q(u)  =  X(j),  <  u  <  j  ,  j  =  1.....N  ,  (1.1) 


where  is  the  jth  sample  order  statistic;  then  given  a  spacing, 

U  =  {u^,...,u^}  (a  set  of  k  real  numbers  satisfying  0<u^<u2<. . .<u^<l) , 

the  ABLUE's  are  linear  combinations  of  the  corresponding  sample  quantiles 

Q(u, ) , . . .  ,Q(u  ).  Explicit  formulae  for  the  coefficients  for  the  Q(u.) 
Ik  i 

may  be  found  in  Ogawa  (1951)  or  Eubank  (1979).  As  these  coefficients 
are  dependent  on  the  form  chosen  for  F,  we  adopt  the  notation  y(U,F), 
ct(U,F)  for  the  ABLUE's  of  y  and  a  in  order  to  indicate  dependence  on 
both  the  distribution  and  the  spacing. 

To  measure  the  performance  of  ah  ABLUE  one  usually  examines  its 
asymptotic  efficiency  relative  to  the  Cramer-Rao  lower  variance  bound. 
Therefore,  assume  that  F  admits  a  density  f  *  F'  and  define  the  quantile 
function  corresponding  to  F  by 

Q(u)  *  inf(x:  F(x)  u}  . 

The  density-quantile  function  is  then  given  by  fQ(u)  =  f(0(u)), 

0  _<  u  _<  1  (c.f.  Parzen  (1979)  for  a  discussion  of  quantile  and  density- 
quantile  functions  and  their  properties) .  Using  this  notation  the 
expressions  derived  by  Ogawa  (1951)  for  the  asymptotic  relative 
efficiencies  (ARE's)  for  the  ABLUE's  may  be  written  as  follows: 
i)  The  ARE  for  simultaneous  estimation  of  u  and  a  is 


ARE(y (U,F) ,  o(U,F)) 


k11(u,f)k22(u,f)-k12(u,f)‘ 


| I (F) |  (1-2) 

where  for  Uq«0,Uj.^=1  and  fQ(0)-fO(l)«*fQ(0)Q(0)-fO(l)Q(l)»0  we  have 
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Kn(U,F) 

K22(u,F) 

k12(u,f) 


k+1  tfQ(ui)  -  fQ(u  )]: 

5  l  ~ 


i*l 


Ui  '  Ui-1 


k+lffQCu^Q^)  -  fQ(ui_1)Q(ui_1)]' 


i-1 


ui  “  Ui-1 


kt1  CfQ(u1)-fQ(ui_1) ] [fQ(ui)Q(ui)-fQ(ut_1)Q(ui_1)] 
i-1  Ui  “  Ui-1 


(1.3) 

(1.4) 

(1.5) 


with 


1(F)  - 
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I  (F) 
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(1.6) 


denoting  the  usual  Fisher  information  matrix  for  location  and 
scale  parameter  estimation  for  the  distribution  F. 

ii)  The  ARE  for  the  estimation  of  u  when  o  is  known  is 

ARE(p(U,F))  -  Ku(U,F)/Iiiy(F)  .  (1.7) 

iii)  The  ARE  for  the  estimation  of  a  when  v  is  known  is 

A 

ARE(a(U,F) )  -  K22(U,F)/Iao(F)  .  (1.8) 

Since  these  ARE’s  are  dependent  upon  the  spacing  for  the  quantiles,  U, 
it  is  common  practice  to  further  optimize  by  selecting  U  to  maximize 
one  of  (1.2),  (1.7)  or  (1.8)  or,  equivalently,  their  numerators 
Ku(U,F)K22(U,F)-K12(U,F)2,  Kn(U,F)  and  K22(U,F)  .  This  is  the  so  called 

optimal  spacing  problem  and  has  been  considered  by  numerous  authors 
(for  references  see  Hassanein  (1977)  and  Eubank  (1981)).  We  wish  to 
consider,  instead,  the  problem  of  selecting  U’s  that  are  robust  in  the 
sense  that  they  provide  high  ARE's  over  several  possible  choices  for  F. 


The  robust  spacing  selection  technique  presented  in  this  paper 
is  motivated  by  an  approach  to  location  parameter  estimation  developed 
by  Chan  and  Rhodin  (1980) .  They  consider  the  problem  of  robustly 
estimating  p  when  F  is  assumed  to  belong  to  a  finite  set  of  probability 
laws,  L.  For  Lei  let  U(L)  denote  the  corresponding  optimal  (maximum 
ARE)  k-element  spacing;  then,  Chan  and  Rhodin  suggest  choosing  a  spacing, 

U(L*),  which  satisfies 

min  ARE(p(U(L*) ,  G))  *  max  min  ARE(p(U(L) ,G)) .  (1.9) 

Ge L  Lei  Gel 

Thus,  one  selects  the  element  of  {U(L);Ltl}  that  is  robust  in  the  sense  of  providing 

the  largest  guaranteed  ARE  (GARE)  of  at  least  min  ARE(p(U  (L*),G))  regardless 

Ge  L 

of  which  law  in  L  generated  the  sample.  In  addition,  this  provides  a 
candidate  for  F,  namely  L*,  and,  hence,  p  may  be  estimated  by  p(U(L*),L*). 

In  this  paper  regression  design  techniques  are  utilized  to  develop 
an  asymptotic  (as  k-*»)  version  of  the  Chan  and  Rhodin  procedure  which, 
in  contrast  to  their  approach,  has  the  advantages  that  it  i)  allows 
for  the  estimation  of  either  or  both  of  p  and  ct  and  ii)  does  not  require 
the  tedious  computation  and  tabulation  of  optimal  spacings.  Due  to 
this  latter  quality,  the  procedure  we  propose  is  amenable  to  use  with 
either  large  or  small  values  of  k.  This  fact  is  of  particular  importance 
since  it  will  usually  be  necessary  to  use  larger  values  of  k  to  insure 
satisfactory  GARE's  than  would  be  required  merely  for  ARE  consideration. 

We  present  our  method  for  spacing  selection  in  Section  2  accompanied 
with  a  discussion  of  its  extension  to  the  problems  of  quantile  estima¬ 
tion  and  parameter  estimation  from  singly  or  doubly  censored  samples. 

A  numerical  example  is  provided  in  Section  3. 


2 .  Robust  spacing  selection 

Parzen  (1979)  has  shown  that,  for  large  N,  location  and  scale 
parameter  estimation  can  be  considered  as  a  regression  analysis  problem 
via  the  model 

fQ(u)Q(u)  «  yfQ(u)  +  ofQ(u)Q(u)  +  OgB(u) ,  ucfO.l]  ,  (2. 

where  aB  *  o/i/W  and  B(*)  is  a  zero  mean  normal  process  with  covariance 
kernel  Cov(B(u^),  B(u2))  *  min  (u^,u2)  -  u^u^  Eubank  (1981)  noted 
that  the  generalized  least  squares  estimators  of  y  and  o  obtained  by 
sampling  from  this  model  at  a  design  U  -  {u^,...,u^}  whose  elements 
satisfy  0<u^<u2<. . .<u^<l  were,  in  fact,  the  ABLUE's.  As  a  result,  it 
was  shown  that  problems  of  spacing  selection  could  be  phrased  as 
regression  design  problems  for  model  (2.1).  The  relationship  between 
these  two  problems  will  be  exploited,  subsequently,  to  select  spacings 
which  have  good  robustness  properties. 

Our  attention  will  be  restricted  to  designs  (spacings)  that  are 
obtained  from  continuous  densities  on  [0,1].  Let  h  be  such  a  density 
with  associated  distribution  and  quantile  functions  H  and  H  ^ .  Then, 
h  generates  a  sequence  of  designs,  {U^K  with 

u  "  jH  » 

i.e.,  {U^}  is  the  design  sequence  whose  kth  element  is  composed  of 
the  (kfl)-tiles  of  H.  For  designs  which  are  obtained  in  this  manner 
it  is  possible  to  characterize  the  asymptotic  (as  k-«°)  behaviour  of 
the  ARE's  of  the  corresponding  ABLUE's  for  the  various  estimation 
situations.  It  follows  directly  from  the  work  of  Sacks  and 
Ylvisaker  (1966,  1968)  and  Eubank  (1981)  that,  under  appropriate 


conditions  on  h: 
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i)  When  both  y  and  o  are  unknown, 

| 1(F) |  -|K(U  ,F) |  =  C  (h,F)  +  o(k-2)  (2.3) 

k  12k  1 

where  K(U^,F)  is  the  matrix  having  ijth  element  K^tU^.F),  i,j*l,2, 
and 

C ,<h,F)  -  /  du  <2.4> 

1  0  h(u) 

with  ip (vi)  defined  by 

<Ku)  =  ( [fQ(u) ] [fQ(u)Q(u) J")1  .  (2.5) 


ii)  When  only  y  is  unknown, 

I  (F)  -  K  n(U.,F)  =  C  (h,F)  +  o(k~2) 

yy  H  K  12k  2  2 

where 

c  (h,F ,  .  /  ilawp!  du  . 

Z  0  h(u)Z 


(2.6) 

(2.7) 


iii)  When  only  a  is  unknown 

Iaa(F)  "  K22(Uk’F)  "  ^2  C3(h»F>  +  °(k‘2> 


where 


c  (h,F) .  /igqwwrv 


du 


(2.8) 


(2.9) 


These  characterizations  can  be  shown  to  hold, for  instance,  when  the 

elements  of  i|>  are  continuous  and  /  (h(u)}  2du  <  00 .  Alternative  conditions 

0 

can  be  found  in  Sacks  and  Ylvisaker  (1968,  Theorem  3.1)  whereas  somewhat 
weaker  restrictions  may  be  deduced  from  the  work  of  Pence  and  Smith  (1981). 
Equations  (2.3),  (2.6)  and  (2.8)  have  the  consequence  that,  for 
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sufficiently  large  k,  what  distinguishes  between  the  efficiency  of  spacings 
selected  from  various  densities  is  their  respective  values  for  the  asymptotic 
constants  C^,  and  C^*  In  particular,  one  can  construct  an  asymptotic 
solution  to  the  optimal  spacing  problem  by  minimizing  C^(h,F)  over  h  for 
a  given  F.  The  solution  was  found  by  Eubank  (1981)  to  be 

(ii/(u)tI  ^(F)i^(u) }^^/X^  ,  when  both  u  and  a  are  unknown, 

2/T 

hp(u)  =*  {[fQ(u)]"}  /X 2  ,  when  only  p  is  unknown  ,  (2.10) 

f[fQ(u)Q(u)]"}  2/3/x  ,  when  only  a  is  unknown, 

3 

where  X^,  X^  and  X^  are  appropriate  normalizing  constants. 

The  elements  of  the  spacing  sequence  generated  by  h  exhibit  the 

r 

same  asymptotic  behaviour  as  a  sequence  of  optimal  spacings  in  the  sense 

that  the  AKE's  for  both  sequences  converge  to  the  same  limit  (either  |l(F)j, 

_2 

I^(F)  or  Iacj(F))  at  the  same  rate  (namely  0(k  )).  In  fact  the  spacings 

obtained  using  (2.10)  have  been  found  to  provide  an  excellent  approximate 

solution  to  the  optimal  spacing  problem  even  for  k  as  small  as  5  or  7 . 

The  reader  is  referred  to  Eubank  (1979,  1981)  for  a  discussion  of  the 

spacings  generated  by  h  and  their  properties. 

r 

The  objective  of  this  paper  is  to  utilize  the  characterizations 

(2.3),  (2.6)  and  (2.8)  to  provide  a  robust  (as  opposed  to  optimal) 

spacing  selection  scheme.  To  do  so,  we  again  utilize  the  constants 

C^(h,F).  However,  instead  of  optimizing  over  h  for  a  fixed  F,  we  will 

optimize  with  respect  to  both  h  and  F  as  follows.  Let  L  denote  a  finite 

set  of  probability  laws  and  assume  that  for  all  combinations  of  L ,GeL 

the  constant  C^Oi^.G)  for  the  estimation  problem  of  interest  is  well 

defined.  An  approximate  (for  large  k)  solution  to  the  problem  of  robust 

spacing  selection  is  then  given  by:  Choose  L *eL  so  that 

max  C. (hT*,G)  ■  min  max  C 
Gel  1  L  LeL  Gei. 


i«V°>  • 


(2.11) 
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r» 


Depending  upon  which  parameters  are  to  be  estimated,  one  then  uses 

either  (u(U*  ,  L*) ,o(U*,L*)) C,  y(U*,L*),  or  a(U*,L*)  for  estimation 
k  k  k  k 

purposes  where  U*  is  the  nth  element  of  the  design  sequence  generated 
k 

by  observe  that  in  the  case  of  location  parameter  estimation 

(i  =  2)  equation  (2.11)  can  be  viewed  as  providing  an  asymptotic  version 
of  (1.9). 

An  illustration  of  the  use  of  (2.11)  is  provided  in  the  next  section. 
First,  however,  we  discuss  some  of  the  computational  details  and  merits 
of  this  approach  as  well  as  certain  of  its  extensions  to  other  robust 
estimation  problems. 

Computation  of  the  constants  (2.4),  (2.7)  and  (2.9)  will  usually 
require  numeric  integration.  For  this  purpose  the  use  of  a  Gaussian 
quadrature  rule  is  recommended.  One  then  constructs  a  table  consisting 
of  the  values  of  C^Or^.G)  for  all  combinations  of  L  and  G  in  L  from 
which  a  law  satisfying  (2.11)  can  be  readily  found.  This  table  has 
the  same  role  in  our  procedure  as  that  of  Tables  3,4,5  and  6  (i.e., 
tabulations  of  the  ARE's  for  optimal  spacings  for  all  L,G  combinations 
corresponding  to  k=2,3,4  and  5  respectively)  in  the  Chan  and  Rhodin 
method.  However,  due  to  the  asymptotic  nature  of  our  solution  a  decision 
reached  using  (2.11)  obtains  for  all  k .  The  computational  savings  derived 
from  this  technique  are,  therefore,  twofold  in  that  i)  a  single  table 
suffices  for  all  values  of  k  and  ii)  the  computation  of  the  constants 
C^,C2  and  C^  requires  only  numeric  integration  as  compared  to  the 
solution  of  nonlinear  equation  systems  required  by  procedures  based 
on  optimal  spacings  (c.f.  Sarhan  and  Greenberg  (1962)  for  examples  of 
such  equations).  Consequently,  criterion  (2.11)  is  amenable  to  use  with 
larger  values  of  k .  The  importance  of  this  fact  should  not  be  over- 
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looked  as,  in  order  to  insure  high  GARE's,  it  will  be  necessary  to  use 
ABLUE's  based  on  more  quantiles  than  would  be  needed  merely  to  obtain 
high  ARE’s. 

Once  an  L*  satisfying  (2.11)  has  been  found  the  next  step  is  to 
compute  the  corresponding  k-point  design,  U£.  Although  the  design 
density  h  ^  will  frequently  not  have  a  closed  form  this  can  still  be 
easily  accomplished  through  numeric  tabulation  of  followed  by 
interpolation  to  find  the  design  points  HL*^k+l^  ’  Tlie 

corresponding  estimator  coefficients  are  then  computed  in  the  usual 
manner  using  the  formulae  found  in  either  Ogawa  (1951)  or  Eubank  (1979) . 

It  should  be  noted  that  whereas  the  spacings  selected  using  (2.11) 
are  robust  this  does  not  imply  that  the  use  of  coefficients 
corresponding  to  L*  is  an  optimal  strategy  for  robustness  considera¬ 
tions  (similar  comments  also  apply  to  the  Chan  and  Rhodin  method). 

In  practice  one  may  wish  to  use  the  spacings  for  L*  in  conjunction 
with  coefficients  for  other  members  of  L  or,  perhaps,  employ  some  type 
of  averaging  procedure.  The  development  of  a  robust  coefficient 
selection  scheme  to  accompany  (2.11)  (as  well  as  (1.9))  is  currently 
an  open  research  problem. 

The  robust  spacing  selection  techniques  presented  here  can  be 

adapted  to  provide  solutions  to  other  robust  estimation  problems. 

For  instance,  it  is  sometimes  difficult  to  interpret  the  parameter 

comparisons  that  are  implied  by  (2.11).  Such  difficulties  can  be 

averted  by  comparing  the  various  location  and  scale  parameter  models 

for  the  laws  in  L  on  the  basis  of  one  or  more  specific  quantiles.  Let 

x  denote  a  specified  percentile  point  and  for  Lei  let  0^  denote  the  corre- 

x— u 

sponding  quantile  function.  The  Tth  quantile  for  the  model  L(— — )  is 


then  u  +  o  Q  (t)  for  which  an  estimator  is  p(U,L)  +  o(U,L)Q  (t). 

L  Li 

The  asymptotic  variance  of  this  estimator  is  found  to  be  proportional 
to  the  trace(tr)  of  K  ^(U,L)  v^djv^fl) t  where  v^d)  =  (1,  Q^(T))t. 

In  a  manner  similar  to  previous  developments  we  have  ,  under 
certain  additional  conditions  on  h  (c.f.  Theorem  4.5  of  Sacks  and 
Ylvisaker  (1968))  ,  that 


with 


tr( [I  A 
C4(h,L) 


(L)-K'1(U,1L)]v(L)v(L)t)  =  -i-j  C.(h,L)  +  o  (n‘2) 
k  T  T  12k*  4 

1  ip(u)tr1(L)vT(L)vT(L)tr1(L)ii-(u) 

=  J  -  -  du  # 

0  h(u) 


(2.12) 

(2.13) 


An  optimal  density  for  this  problem  is  proportional  to 
4(«)tI~1(L)v  (L) v  (L)tr1(L)»|)(u)}1/3 

T  T 

and,  hence,  (2.13)  can  now  be  used  in  (2.11)  to  determine  robust 
spacings  for  the  purpose  of  quantile  estimation. 

It  is  also  possible  to  modify  (2.11)  to  provide  robust  spacings 
for  problems  of  estimation  from  left  and/or  right  censored  samples 
such  as  those  considered  by  Cheng  (1980) .  This  can  be  accomnlished 
through  use  of  the  optimal  densities  for  spacing  selection  from  censored 
samples  given  in  Eubank  (1981). 


3.  Numerical  example. 

As  an  illustration  of  the  use  of  the  spacing  selection  technique 
provided  in  Section  2  consider  the  data  consisting  of  lifetimes  for 
417  40-watt  internally  frosted  incandescent  lamps  presented  in 
Davis  (1952) .  This  data  has  also  been  analyzed  by  Chan  and  Rhodin  (1980) 
under  the  assumption  that  the  parent  distribution  was  either  a  member 
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of  (or  well  approximated  by)  the  Tukey  lambda  family  of  distributions 
or  was  a  normal,  double  exponential  or  Cauchy  distribution.  Through 
goodness  of  fit  considerations  they  chose  L  to  consist  of  the  normal 
distribution  and  those  members  of  the  Tukey  lambda  family  having  shape 
parameters  -.1,  0,  .1  and  .14.  To  facilitate  comparison  with  their 
results  L  will  be  chosen  similarly  here.  Letting  L(X)  denote  the 
Tukey  lanbda  distribution  with  shape  parameter  X,  we  take  1={L(-.1) ,L(0) , 
L(.l),  L(.14)}  where,  for  simplicity,  the  normal  has  been  omitted  since 
it  is  approximated  by  L(.14). 

The  quantile  function  for  L(X)  is 
Qx(u)  -  X_1[uX  +  (l-u)XJ 
with  associated  density-quantile  function 
(fq)x(u)  =  [ux_1  +  (i-u)x_1]_1. 


Thus,  the  optimal  density,  =  ^X*  ^or  locatic‘  parameter  esti¬ 

mation  is  proportional  to 


{(fQx)"(u)}2/3=(uX_1+(l-u)X_1)'2{(X-l)X[uX  2-(1-u)X'2]2-(X-1)(X-2)[u(1-u)]X  3'2/3. 


As  ^^(u)  =  h^Cl-u),  the  spacings  generated  by  h^  will  be  symmetric,  i.e., 
u^  *  l-u^+^  This  has  the  consequence  that,  u(Uk,  L(X))  is  scale 
invariant  (c.f.  also  Chan  and  Rhodin  (1980)  pg.  228  for  further  comments 
on  this  property) . 

The  constants  C2(hx>  L(n))  =  C2(X,n)  have  the  form 
C2(X,n)  *  /1[(fQn)"(u)/hx(u)]2du. 

The  values  of  C2(X,n)  have  been  tabulated  for  all  X,n  e  (-.1,0, .1, .14} 
and  are  provided  in  Table  1.  It  follows,  upon  examination  of  Table  1, 
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that  the  minimum  value  of  max  C9(X,n)  occurs  at  X*  =  .14  and,  hence, 
h  ^  generates  spacings  that  are  robust  relative  to  L.  It  is  of  interest 
to  note  that  this  choice  of  X  agrees  with  the  one  made  by  Chan  and 
Rhodin  for  k=2  and  is  their  second  choice  (i.e.  has  the  second  largest 
minimum  ARE)  when  k=3.  It  is  unfortunate  that  tables  are  not  available 
for  larger  values  of  k,  such  as  k -  7,9  or  perhaps  even  20,  where  one 
might  expect  the  results  from  these  two  techniques  to  consistently 
coincide . 

Numeric  tabulation  of  the  density  h  ^  reveals  that  the  robust  5 
element  spacing  is  U,-  =  {.0105,  .1628,  .5,  .8372,  .9895}  from  which 
the  corresponding  estimator  is  computed  to  be  (using  Ogawa's  formula) 

u(U5,L(.14))  =  .0396[Q( ,0105)+Q( .9895) ]  +  . 2585[Q( . 1628)-K}( . 8372) ] 

+  .4037Q( . 5) 

=  .0396[609+1550]  +  . 2585 [883  +  1225]  +  .4037[1037] 

=  1049.0513. 

This  is  to  be  compared  with  the  value  1044.3  obtained  by  Chan  and 
Rhodin  (1980)  also  through  the  use  of  a  5-quantile  ABLUE. 

As  larger  values  of  k  will  usually  be  required  to  insure  suffi¬ 
ciently  high  GARE's  it  is  important  to  note  that  once  h  ,,  has  been 

.  14 

tabulated  spacings  of  any  desired  size  can  be  easily  computed.  For 
instance,  when  k  =  7  the  robust  spacing  is  *  {.0034,  .0574,  .2353, 

.5,  .7647,  .9426,  .9966}  with  corresponding  coefficients  {.0148,  .0995, 
.2374,  .2967,  .2374,  .0995,  .0148}.  The  resulting  estimator  is  then 
found  to  be  u(U?,L( .14))  =  1045.97. 
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It  is  also  of  interest  to  ascertain  how  U,.  and  U^,  perform  in 
terms  of  GARE.  Consequently,  the  ARE’s  for  these  spacings  have 
been  computed  over  the  various  laws  in  L  and  are  provided  in  Table  2. 
The  minimum  or  guaranteed  ARE  when  k  =  5  is  found  to  be  . 908  whereas 
when  k  =  7  the  GARE  is  .9466.  Although  the  set  L  considered  by  Chan 
and  Rhodin  includes  the  normal  while  ours  does  not,  it  is  still 
noteworthy  that  the  GARE  for  is  quite  similar  to  the  GARE  of 
.9015  for  the  robust  5-quantile  ABLUE  based  on  optimal  spacings 
utilized  by  Chan  and  Rhodin. 
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TABLE  2 


ARE's  for  and  for  the  various  laws  in  L  . 


n 

ARE(u(U5>L(n))) 

ARE(u(U7,L(tO)) 

-.1 

.9080 

.9466 

0 

.9163 

.9513 

.1 

.9156 

.9506 

.14 

.9111 

.9474 

I 
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